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ABSTRACT 

We derive a simple approximate model describing the early, hours to days, UV/optical supernova 
emission, which is produced by the expansion of the outer < 1O -2 M0 part of the shock-heated enve- 
lope, and precedes optical emission driven by radioactive decay. Our model includes an approximate 
description of the time dependence of the opacity (due mainly to recombination), and of the deviation 
of the emitted spectrum from a black body spectrum. We show that the characteristics of the early 
UV/O emission constrain the radius of the progenitor star, its envelope composition, and the 
ratio of the ejecta energy to its mass, E/M. For He envelopes, neglecting the effect of recombina- 
tion may lead to an over estimate of i?* by more than an order of magnitude. We also show that 
the relative extinction at different wavelengths (A\ — Ay) may be inferred from the light-curves at 
these wave-lengths, removing the uncertainty in the estimate of i?» due to reddening (but not the 
uncertainty in E/M due to uncertainty in absolute extinction). The early UV/O observations of the 
type lb SN2008D and of the type Hp SNLS-04D2dc are consistent with our model predictions. For 
SN2008D we find i?* w 10 11 cm, and an indication that the He envelope contains a significant C/O 
fraction. 

Subject headings: shock waves — radiative transfer — relativity — supernovae: general — supernovae: 
individual (SN 2008D) — stars: evolution 



1. INTRODUCTION 

During the past four years, the wide field X-ray de- 
tectors on board the Swift satellite enabled us, for the 
first time, to "catch" supernova (SN) explosions very 
close to the onset of the explosion. In two cases, 
SN2006aj and SN2008D, the usual optical SN light curve 
was observed to be preceded by a luminous X-ray out- 
burst, which has triggered the X-ray detectors, fol- 
lowed by a longer, ~ 1 day, duration UV/ O emission 
(|Campana et al.ll2006b ISoderberg et"aT1 120081 ) . Analysis 
of the later optical SN emission revealed that both were 
of type Ib/c, probably produced by compact Wolf-Rayet 
(hereafter WR) pro genitor stars, which lost most of their 
Hydrogen env e lope (lPian"eT al. 2006; Maz zali et alj2006t 
Modiaz etH\ 120061; IMaeda et all 120071: IMazzali et all 



20081: iMalesani et al.N2009| ) Following the detection of 



the early UV/O emission from these SNe, a search was 
conducted for UV/O emission from SNe that fall within 
the deep imaging survey of the Galaxy Evolution Ex- 
plorer (GALEX) space telescope, leading to the detec- 
tion of ea rly (~ lday) rising UV/O emission for two type 
II-p SNe (jGezari et al]l200a ISchawinski et al.l[2008h . for 
which the progenitors are likely red super giants (here- 
after RSG) with large Hydrogen envelopes. 

X-ray outbursts followed by early UV emission have 
long been expected to mark the onset of SN explosions. 
The SN shock wave, which travels through, and ejects, 
the stellar envelope, becomes radiation-mediated when 
propagating through the envelope (for review see, e.g., 
iWooslev fc Weaverlll986l ). As the shock propagates out- 
ward, the (Thomson) optical depth of the plasma ly- 
ing ahead of it decreases. When this optical depth be- 
comes comparable to the shock transition optical depth, 
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t s ~ c/v s (y s is the shock velocity), the radiation es- 
capes ahead of the shock. T his leads to an expected 
"shock breakout" X-ray flash (|Colgatd 11971 lFalldll978t 
Klein fe Chevalierlll978l) lasting for 10's to 100's of sec- 
onds. Following breakout, the stellar envelope expands 
and cools (nearly adiabatically) . As the photosphere 
penetrates into the outer shells of the envelope, the (adia- 
batically cooled) radiation stored within the envelo pe es- 
capes , leading to an expected early UV/O emission 
1978). In this paper we focus on the early , ~ 1 day, part 
of this UV/O emission (although it may dominate the to- 
tal emission for much longer, e.g. in type II-P SNe). 

The interpretation of SN associated X-ray outbursts 
as due to "shock breakout" is not generally accepted, 
due mainly to the fact that while a ~ 0.1 keV 
thermal spectrum was expected, the observed spectra 
are non thermal and extend b eyond 10 keV. Some 
authors dCampana etaLl 120061: IWaxman et al.1 120071 : 
ISoderberg et al.l 120081 ) argued that this is due to shock 
breakout physics. Others argued that the X-ray bursts 
can not be explained within this framework and im- 
ply the existence of relativistic energet ic jets penetrat- 
ing through the stellar mantle /envelope (ISoderberg et all 
20061: iF^jTeTall [2lM iGhisellini et all I2007t ILil l2007t 
Mazzali et al.l 120081 : Hi 120081 ). A recent derivation of 



the structure of mi ldly and highly re lativistic radiation- 
mediated shocks (jKatz et al.l [2010) shows that fast, 
v s /c > 0.2, radiation mediated shocks produce photons 
of energy far exceeding the ~ 0.1 keV downstream tem- 
perature, reaching 10's to 100's of keV. This suggests 
that the observed outbursts may indeed be due to shock 
breakout (Note, that for SN2006aj there is an additional 
challenge, beyond the X-ray spectrum: The energy in- 
ferred to be deposited in the fastest part of the ejecta 
far exceeds that expected from shock acceleration in the 
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envelope) . 

Our focus in the current paper is not on the X-ray 
outburst, but rather on the early, ~ 1 day, UV/O emis- 
sion that follows it. Model predictions for the early 
UV/O emission were derived mainly using numerical 
calculations (e.g. iFalkl Il978t IE nsman fc Burrows! [1992; 
iBlinnikov et~afll2000t l2002t iGezari et al.ll2008l )7~Follow- 
ing the detection of SN2006aj, an an alytic model has 
been constructed l|Waxman et al.ll2007t) . One of the main 
advantages of the analytic model is that it provides ex- 
plicit analytic expressions for the dependence of the emis- 
sion on model parameters, thus making both the use of 
observations for determining parameters and the iden- 
tification of model uncertainties much easier and more 
straightforward. It was shown, in particular, that the 
photospheric temperature of the expanding envelope de- 
pends mainly on the progenitor's radius and on the opac- 
ity, T p h approximately proportional to rV^ , and that the 
luminosity L is approximately proportional to [E /M)R*. 
This implies that the progenitor radius, which is only 
poorly constrained by other observations, and E/M may 
be directly determined by measuring and analyzing the 
early UV/O emission. 

In should be noted here that most numerical models, 
which are used for analyzing SN light curves, focus on 
the long term radioactively driven emission. Such mod- 
els do not describe the early UV/O emission (e.g. figure 
3 in lTanaka et al. 2009), largely due to the fact that they 
lack the resolution required to properly describe the evo- 
lution of the outer ~ 10~ 2 M Q part of t he shock-heated 
envel ope, which drives the early emission (Waxman et al. 
[2007J) ■ Numerical calculations that allow a proper treat- 
ment of the early emission are, on th e other hand, com- 
putationally very dem anding (see e.g. IGezari et al.ll2008l 
IBlinnikov et al.ll2000D . It is difficult to use these mod- 
els for obtaining the dependence of predictions on model 
parameters and therefore for using observations to con- 
strain these parameters. 

In the near future, we expect an increasing rate of de- 
tection of early UV/O SN emission. Ground based SN 
sur veys with high rate sampling (le ss or order of a day, 
e.g. lLaw et al.l 120091: lOuimbvi I2006D will detect SNe at 
the early stages of their expansion, with a bias towards 
detection of the more abundant type II SNe. The de- 
tection of X-ray outbursts, that were observed to mark 
the onset of several SN Ib/c explosions, suggest that the 
early UV/O emission from these SNe may be detectable 
by follow-ups of X-ray triggers. The MAXI (Monitor 
of All-sk y X-ray Image) exper iment on board the Kibo 
module ()Matsuoka et al.lll997l ). which was launched this 
year, is expected to detect events similar to SN2008D 
at a rate of up to a few per year. The EXIST satellite 
(jGrindlav et al.ll2003t iBand et al.M2008D . which is still in 
planning, will further increase the event rate and improve 
the X-ray spectral coverage (the detection rate may also 
be enh anced by smaller dedicated X-ray observatories, 
see e.g. [Calzavara & Matzner 2 0041) . 

The main go al of the current pap er is to extend the 
analytic model (Waxm an et al.| [2007) to include a more 
realistic description of the opacity and its variation with 
time (mainly due to recombination), and to include an 
approximate description of the deviation of the emitted 
spectrum from a black body spectrum (due to photon 



diffusion). Our results will facilitate the use of upcoming 
observations of the early emission from SNe for constrain- 
ing progenitor and explosion parameters. We consider 
several types of progenitor envelopes: Dominated by Hy- 
drogen, as appropriate for RSG and BSG progenitors, as 
well as envelopes dominated by He, C/O, and He-C/O 
mixtures representing different de grees of H/He strip- 
ping, due to wind mass l oss (e.g. iWooslev et al.l [T993D 
or binary interaction (e.g. iNomoto et al.lll995() . and dif- 
ferent e volution scenarios (e.g. due to rotation induced 
mixing, iMevnet fc Maederll200l iCrowtherlfeOOTf ). 

For completenes s, we first present in § [2] the simple 
model derived in ( Wa xman et al.l [20071) . The main as- 
sumptions adopted are that the envelope density drops 
near the stellar edge as a power of the distance from the 
edge, eq. ((TJ), tha t the SN shock velocity m ay be approxi- 
mated (following IMatzner fc McK cc 1999) by an interpo- 
lation between the Sedov-von Neumann- Taylor and the 
Gandcl'Man-Frank-Kamenetskii-Sakurai self-similar so- 
lutions, eq. ([2]), and that the opac ity, k , is space and time 
independent. We also show, in § 12.31 that the effects of 
p hoton diffusion on the predicted luminosity ( considered 
in lChevalierlll992tlChevalier fc Frans son 2008) are small. 
The model is extended in § [3] to include a more realistic 
description of the opacity and an approximate descrip- 
tion of the effect of photon diffusion on the spectrum. In 
§[4] we show that the relative extinction at different wave- 
lengths (A\ — Ay) may be inferred from the light-curves 
at these wave-lengths. In §[5] we compare our model pre- 
dictions to observations of the early emission available for 
two SNe, arising from RSG and BSG progenitors, and to 
detailed numerical simulations that were constructed to 
reproduce these observations. In §[6] we use our model to 
analyze the early UV/O observations of SN2008D. Our 
main results are summarized and discussed in § [71 

2. A SIMPLE MODEL 

We first derive in § 12.11 the density, velocity and tem- 
perature profiles of the (post-breakout) expanding stellar 
envelope. We then derive in § !2.2l the radius and temper- 
ature of the photosphere. These derivations are carried 
out under the simplifying assumption, that photon dif- 
fusion is negligible (and hence that the flow is adiabatic) 
below the photosphere. This assumption is justified in 



2.1. Expanding ejecta profiles 

The UV/O emission on a day time scale arises from 
the o uter < 10 _2 M Q shell of the ejecta (Waxman et al. 
I2007h . Neglecting the shell's self-gravity and its thickness 
(relative to -R*), the pre-explosion density profile within 
the shell may be approximated by (|Chandrasekharfl939f ) 



Po(m) = Pi/2<5" 



(1) 



where S = (1 — ro/i?*), r$ is the radius, and 
n = 3,3/2 for radiati ve and efficiently c o nvecti ve en- 
velopes respectively. Matzner & McKee (1999) have 
shown that the velocity of the SN shock within 
the envelope is well approximated by an interpola- 
tion between the Sedov-von Neumann-Taylor and the 
Gandel' Man-Frank-Kamenetskii-Sakurai self similar so- 
lutions, (|Von Neumann! 119471: ISedovl 119591 : iTavloii H950I : 



IGandel'Man fc Frank-Kamenetskii|[T95l ISakurailll9"60h 
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where E is the energy deposited in the ejecta, m(ro) 
is the ejecta mass enclosed within ro, A v ~ 0.79 and 
fa ~ 0.19. For r -> i?*, eq. fl5) reduces to 
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where / p = pxi^j 'p and p is the average ejecta density, 
p = 3M/4nRl. iCalzavara fc Matznerl (|2004 Appendix 
A) have derived the values of f p expected for various pro- 
genitors (using their notation, f p — (3/47r)pi/p*). For 
BSGs they find f p varying nearly linearly with mass, 
from 0.031 to 0.062 for 8.5 < M/M Q < 18.5, and for 
RSGs they find 0.079 < f p < 0.13. We show below that 
the results are not very sensitive to the value of f p . 

In what follows we replace the Lagrangian coordinate 
ro with 



6 m (6) = M- 



dr4:Trr 2 po(r) 



3/ P 



rn+l 



(1-8) R. 



1 



(4) 



the fraction of the ejecta mass lying initially above ro. 

As the radiation mediated shock passes through a fluid 
element lying at ro, it increases its pressure to 



6 2 
Po = ljP0V s 



(5) 



and its density to 7po (recall that the post shock en- 
ergy density is dominated by radiation). As the shocked 
fluid expands, it acceler ates, converting its intern al en- 
ergy to kinetic energy. iMatzner fc McKed (|1999ft have 
shown that the final velocity, Vf(ro), of the fluid initially 
lying at ro is well approximated by Vf(ro) = fv(fo)v a (ro) 
with /„ « 2. The value of f v depends on the curvature 
of the shells. The effect of this dependence is considered 
below. 

Equations ([3]) and (O hold as long as the shock width 
is much smaller than the width of the stellar envelope 
shell lying ahead of the shock. When the shock reaches 
a radius at which the optical depth of the shock tran- 
sition layer, t s = c/v s , becomes comparable to the 
optical depth of the shell lying ahead of the shock, 
r i£ MSmK/AnRl, the radiation "escapes" ahead of 
the shock, producing a "shock breakout flash", and the 
shock can no longer be sustained by radiation. The 
mass fraction at which breakout takes place is (e.g. 
IMatzner fc McKed fl999h 
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2 x 10" 



f -0.07^2.3 
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Elf{M/M Q ) 



0.57„1.1 
^0.34 



(6) 



For smaller values of 5 m , S m < (5 m ,BO> the velocity and 
pressure profiles are shallower than given by eqs. ([3]) 
and ©. 

Assuming adiabatic expansion, neglecting photon dif- 
fusion below the photosphere, the pressure and density 
of the expanding fluid are related by 



p{$ m ,t) 



p(S m ,t) 



7po(<5m) 



4/3 



Po(8 m ). 



(7) 



Once a fluid shell expands to a radius significantly larger 
than , its pressure drops well below p and its velocity 
approaches the final velocity Vf. At this stage, Vft^R*, 
the shell's radius is approximately given by 



r(S m ,t) = v f (5 m )t 
and its density is then given by 



M 



A-rtrH V dS, 



dv 



1 M 



fa n 4nt 3 v 3 f 
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(9) 



The resulting density profile is steep, dlnp/dlnr = 
dlnp/dlnvf = -3 - (n + I)/ fan » -10. 

2.2. Photo spheric temperature and radius 

For a time and space independent opacity n (which 
applies, e.g., for opacity dominated by Thomson scatter- 
ing with constant ionization), the optical depth of the 
plasma lying above the shell marked by S m is 

• dSL 



Jr(5 m ,t) 4 7T JO r {°m) 

1 kM6„ 



1 



2fan/(l + n) Airt 2 v}(6 m )' 



(10) 



where the last equality holds when eq. ([8]) is satisfied. 
We define the Lagrangian location of the photosphere, 
S m , P h, by T(S m = 8 m , p h,t) = 1. We consider two type 
of envelopes: radiative envelopes typical to blue super 
giants (BSG), for which we take n = 3, and efficiently 
convective envelopes typical to red super giants (RSG), 
for which we take n = 3/2 (note, that inefficient con- 
vection may lead to a more complicated density profile). 
For ft = 3/2 and n = 3 envelopes we have 



170.81 o 

8 m , ph (t) = 2.4 x 10-3 / -o.i2 ^^i_ ,i.63 ( „ = }) 



6 m , P h(t) = 2-6 x 10- 3 /p" ' 073 



(M/M )i-««§ : g 



(m/m )i-««8:S 



t^ (» = 3), 
(11) 



10 5 i5 s, and k 



where E — 10 5 E 51 erg, t = iu is 
0.34kq 34 c m 2 / g. Here, and in w hat follows, we use 
( following IMatzner fc McKed" [1991 fa = 0.1909, /„ = 
2.1649, and A v = 0.7921 for n = 3/2 and fa = 
0.1858, /„ = 2.0351, and A v = 0.8046 for n = 3. Us- 
ing eqs. (J3J, ©, and {?]), we find that the radius and the 
effective temperature of the photosphere are given by 
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cm (n 



r ph (i)=3.3xl0 14 / p " 
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51 °- 34 '°- 78 cm(n = 3), 
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and 



T ph (i) = 1.6/. 
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p0.016 p 1 / 4 
^51 ■ ft *,13 
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eV (n = 3). 
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Here, i?* = 10 13 i?*,i3 cm. The dependence on n and 
on f p is weak. Note, that eq. (fT3|) corr ects a typo (in 
the n umerical coefficient) in eq. (19) of (iWaxman et al.l 
[2003). 

As mentioned in the introduction, the photosphcric 
temperature is weakly dependent on E and M and is 
approximately linear in (R^/k) 1 / 4 . The photospheric ra- 
dius, on the other hand, does not depend on i?„ , is weakly 
dependent on k and is approximately linear in E 0A /M 03 . 
The luminosity predicted by the simple model described 
here, L = 47rcrr 2 , h T p 4 h , is 



8.5 x 10 4 



p0.92 p 



/O^(M/M )O-84 K ; 



0.84^0.92^5 



-0.16 



0.34 



for n — 3/2, and 



L = 9.9 x 10 



42 



P0.85n 
£j 51 SX 



*,13 



-0.35 



ergs 1 (14) 



ergs" 1 (15) 



for n = 3. 

Our simple description of r p h T p h, and L, eqs. (fl"2UT5|) . 
holds for 5 m p h > i m BO. i- e - f° r (comparing eqs. © 
and jnj) 



i > £ B o = 0.05 x 10 



5 (M/Mq)°- 6 R} a 



*,13 



0.2 
^0.34 



P0.9 
^51 



(16) 



This requirement also ensures that the ejecta shells have 
expanded and cooled significantly, and thus reached their 
terminal velocity Vf. For these times, the approximation 
of eq. (|7|) holds, and the radius of each shell is well ap- 
proximated by r = t Vf . 

The value of /„ deviates from 2 for large 6 m , due to the 
incre asing curvature of the shells (e.g. Matzne r fc McKeel 
1999). Requiring the deviation not to exceed 30%, in 
which case the error in eq. (| 13[) is smaller than 15%, 
implies limiting the analysis to times 

t < (f p /0.07)°- 69 (M/M e )K^ M E^- 5 x t lv , (17) 

where t fv = 4.5 x 10 6 s, 1.2 x 10 5 s for n = 3/2 and n 
,'i envelopes respectively. Thus, the curvature effect is 
negligible on day- week time scale for M ~ 10M©, and 
may be significant on day time scale only for low mass 
ejecta. 

2.3. Photon diffusion 

Let us next examine the assumption, that photon dif- 
fusion does not lead to strong deviations from adiabatic 
expansion below the photosphere. The size of a region 
around r(S m , t) over which the diffusion has a significant 

effect is D(8 m ,t) ~ y/ ct/3kp(5 m , t). Thus, the radius 
I'd = f(8m,di t) above which diffusion affects the flow sig- 
nificantly may be estimated as D(S m — 5 m>c t, t) — r(S m = 
S m ,d,t). This gives 



il4 £ — 0.069 



E, 



0.45^.0.1 
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^0.34 j.0.79 



r d (t) = 3.7 x 10"/, (M/Ms)0 , 5 - 5 



f>- 79 cm(n = ^), 



r d (t) = 3.8 x 10 14 / P " 



0.44^0.12 



-^51 K > 



(M/Mq) 



°ff^- 75 cm(n = 3). 



(18) 

This radius is similar, and somewhat larger than, the 
photospheric radius given by Eq. (TT21) . The rapid in- 



crease of the diffusion time, ~ 3npr 2 /c, at smaller radii 
implies that diffusion does not significantly affect the 
fluid energy density below the photosphere. Next, we 
note that in regions where the diffusion time is short, 
the luminosity carried by radiation, L oc r 2 dp/d,T, is ex- 
pected to be independent of radius. The steep depen- 
dence of the density on radius (<iln p/dhir ~ —10) then 
implies that the energy density roughly follows pocr, i.e. 
T oc r 1 / 4 . This temperature profile is close to the adia- 
batic profiles derived in § 12.11 for which T oc t °- 28 ,t - 27 
for n — 3, 3/2 respectively. Thus, diffusion does not lead 
to a significant modification of the pressure and temper- 
ature profiles also at radii where the diffusion time is 
short. 

The validity of the above conclusio ns may be teste d 
by using the self-similar solutions of iChevalierl (| 1992f> . 
which describe the diffusion of photons in an expand- 
ing envelope with a density following p oc r - m t m -3 
and initial pressure p oc r t l ~ A . The evolution of the 
ejecta density and pressure derived in § 12.11 follows, for 
v f t > p = Br- m t m ~ 3 with m - 3 = (1 + n)/nPu 
and p = Ar~ l t l ~ 4 with I = (37 — 2) + (7 + n) / n / 3i , wh ere 
7 = 4/3. Applying the solutions of IChevalierl fl992) to 
these profiles one finds the same rj, as given by eq. (fT5j) 
and 
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4irr 2 h o-T 4 , differ from those derived in 



. (. 2 °) 

for n = 3/2 (The parameter q of the self-similar 
solutions is q = 0.495,0.683 for n = 3 and n = 3/2 re- 
spectively; The density and pressure coefficients are A = 
2.95E 3 - 89 R*f- 0A4 M- 2 - 89 ,53.7E 4 - 95 R*f-°- 89 M- :i - 95 
and B = i00E 3 - 59 f- a - 33 M- 2 - 59 ,10 3 E 4 - 37 f- a e7 M- 3 - 37 
for n = 3 and n = 3/2 respectively). 

The parameter dependence of L c derived using the self- 
similar diffusion solutions is similar to that obtained by 
the simple model of § 12.21 and the normalization of L c 
derived using the self-similar diffusion solutions differs 
from the results of §H eqs. (jTSJ and (JTU), by w 10%. 
The effective temperatures derived from the diffusion so 
lutions, via L 

the previous section by 1% to 5% 

Thus, the effects of diffusion on the luminosity and 
on the effective temperature are small, as expected. It is 
important to emphasize in this context that since the dif- 
fusion approximation breaks down near the photosphere, 
the results obtained using the self-similar diffusion solu- 
tions are not necessarily more accurate than those de- 
rived (e.g. in the previous section) by neglecting photon 
diffusion below the photosphere. Improving the accuracy 
of the simple model requires a transport, rather than a 
diffusion, description of the photon propagation. The dif- 
ferences between the results obtained neglecting diffusion 
and including it may be considered as a rough estimate 
of the inaccuracy of the model. 

Finally, th e follow ing note is in place here. The results 
of IChevalierl (|1992D for L and r d (eqs. 3.19 and 3.20) 
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are different both in normalization and in scaling from 
those derived here. This is due to some typographical 
er rors in earlier eqs. of that paper. When corrected, 
in lChevalier fc Franssonl (|2008l ). the results obtained us- 
in g the diffusion s o lution s are similar to those obtained 



Waxman et al.l J2007) and here. The difference in 



the numerical co efficient of the photos pheric tempera- 
ture, eq. (19) oflWaxman et all (|2007D and eq. (5) of 
iChevalier fc Fransson (120081) . is due mainly to the typo 
in eq. (19) of iWaxman et all (|2007| ). which is corrected 
in eq. (fl3|) above. 



3. MODEL EXTENSIONS 

The approximation of space and time independent 
opacity is justified at early times, when the envelope is 
highly ionized and the opacity is dominated by Thom- 
son scattering. On a day time scale, the temperature of 
the expanding envelope drops to ~ 1 eV, see eq. (fl3|) . 
At this temperature, significant recombination may take 
place, especially for He dominated envelopes, leading to 
a significant modification of the opacity. The model pre- 
sented in § [2] is generalized in § 13.11 to include a more 
realistic description of the opacity. The deviation of the 
emitted spectrum from a black body spectrum, due to 
photon diffusion, is discussed in § 13.21 A brief discussion 
of the effect of line opacity enhancement due to velocity 
gradients is given in § 13.31 Throughout this section, we 
use the density structure given by eq. (fTJ with n = 3, 
as appropriate for radiative envelopes. As explained in 
the previous section, the results are not sensitive to the 
exact value of n. 

3.1. Varying opacity 

In order to obtain a more accurate description of the 
early UV/O emission, we us e the mean op acity provided 



eq 



ly UV/U emission, we use the mean opacity provided 
the OP project tables (jSeatonl 120051 ). We replace 
HOI) with 



r(S m ,t) 



dr P K[T(5 m ,t),p(S m ,t)], (21) 



r(<W) 



where k(T, p) is the Rosseland mean of the opacity, and 
solve r{8 m — 8 m ^ p h, t) = 1 numerically for the location 
of the photosphere. In order to simplify the comparisons 
with the suggested analytical models, in the reminder 
of this section we shall take the ejecta properties in the 
limit of eq. (0. 

3.1.1. H envelopes 

Consider first explosions in H dominated envelopes. In 
figure [1] we compare the temperature of the photosphere 
calculated using the OP tables with the results given by 
eq. (|13|) for k = 0.34cm 2 /g, corresponding to fully ionized 
H. The difference in T p h obtained by the two methods 
is smaller than 10% for T p h > 1 eV. At lower temper- 
atures, the k = 0.34cm 2 /g approximation leads to an 
underestimate of Tph, by « 20% at 0.7 eV. This is due to 
the reduction in opacity accompanying H recombination. 
The reduced opacity implies that the photosphere pen- 
etrates deeper into the expanding envelope, to a region 
of higher temperature. The photospheric radius is not 
significantly affected and is well described by eq. (fl2)) . 
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Fig. 1. — Photospheric temperature (top panel) and the ratio of 
color to effective (=photospheric) temperatures (bottom panel) for 
explosions in H dominated envelopes (H with solar composition of 
heavier elements). The top panel compares the analytic approxi- 
mation of eq. 11 31 1 for fixed, k = 0.34cm 2 /g, opacity (dashed lines) 
with the numerical solution of eq. 12111 using OP table opacities 
(solid lines). The color to effective temperature ratio was calcu- 
lated as explained in § 13.21 Results are shown for E = 10 51 erg, 
M = 1 Mq, and three progenitor radii, {30, 10, 3} X 10 12 cm. 



3.1.2. He envelopes 

Let us consider next explosions in He dominated en- 
velopes. In this case, the constant opacity approximation 
does not provide an accurate approximation for T p h- We 
therefore replace eqs. ([12")) and ([13")) with an approxima- 
tion which takes into account the reduction of the opac- 
ity due to recombination. On the time scale of interest, 
hour < t < day, the photospheric temperature is in the 
energy range of 3eV > T > leV. In this temperature 
range (and for the characteristic densities of the photo- 
sphere), the opacity may be crudely approximated by a 



6 



broken power law, 
K = 0.085 cm 2 g 



_! f (T/1.07 eV) 88 
\ (T/1.07 eV) 10 



T > 1.07 eV 
T < 1.07 eV 



(22) 



Using this opacity approximation, we find that cq. ([T 
for the photospheric temperature is modified to 

ri.33cv/-°-°x?st 5 - 

phW \1.07 eV(t/ti)-°- ia 



0.38 



T ph > 1.07 eV 
T ph < 1.07 eV ' 

(23) 

Here, i?* = 10 12 i?» j i2 cm and if, is the time at which 
r p h = 1.07 eV, and we have neglected the dependence 
on E and M, which is very weak. The photospheric 
radius, which is less sensitive to the opacity modification, 
is approximately given by 



-0-038^0.39 (M/Mq) - 



°- 28 ^ 75 cm. 



r ph (t) = 2.8 x 10 14 /; 

Here we have neglected the dependence on R*, which is 
weak. For T p h > 1.07 eV, the bolometric luminosity is 
given by 



L = 3.3 x 10 42 



/o-is(M/M )o 



07 



t 5 -°- 03 ergs- 1 . 



(25) 



The deviation of /„ from the assumed value of 2 for large 
6 m (discussed in § I2.2p has only a negligible effect on 
eqs. (J23J) and ([25]). 

In figure [2] we compare the approximation of eq. (|23[) 
for Tph with a numerical calculation using the OP opac- 
ity tables. The approximation of eq. (|2"5|) holds to better 
than 8% down to Tph — 1 eV. The temperature does not 
decrease significantly below ~ 1 eV due to the rapid de- 
crease in opacity below this temperature, which is caused 
by the nearly complete recombination. 

The following comment is in place here. The strong 
reduction in opacity due to He recombination implies 
that the photosphere reaches deeper into the envelope, 
to larger values of 5 m . The plus signs in fig. [2] denote 
the time at which S miP h = 0.1. For such a large mass 
fraction, the initial density profile is no longer described 
by eq. ([T]) and the evolution of the ejecta is no longer 
given by the eqs. of § 12.11 This further complicates the 
model for the emission on these time scales (see § 16.11 for 
further discussion). 

3.1.3. He-C/O envelopes 

Let us consider next envelopes composed of a mixture 
of He and C/O. At the relevant temperature and density 
ranges, the C/O opacity is dominated by Thomson scat- 
tering of free electrons provided by these atoms, and is 
not very sensitive to the C:0 ratio. Denoting by 1-Z the 
He mass fraction, the C/O contribution to the opacity 
may be crudely approximated, within the relevant tem- 
perature and density ranges, by 



k = 0.043 Zcm 2 g^ 1 (T/l eV) 



1.27 



(26) 



This approximation holds for a 1:1 C:0 ratio. How- 
ever, since the opacity is not strongly dependent on this 
ratio, Tph obtained using eq. (f2l)]) (eq.[57]) holds for a wide 
range of C:0 ratios (see discussion at the end of this sub- 
section). At the regime where the opacity is dominated 
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Fig. 2. — Photospheric temperature (top panel) and the ratio 
of color to effective (=photospheric) temperatures (bottom panel) 
for explosions in He dominated envelopes (Heli um mass fraction of 
98%, C/N/O/Ne fractions taken from fig. 18 of lMevnet fc Maederl 
2003^ The top panel compares the analytic approximatio n of 
eq. d23l ) (dashed lines) with the numerical solution of cq. [12 1 t i 
using OP table opacities (solid lines). The colo r to effective 
temperature ratio was calculated as explained in § 13.21 Results 
are shown for E = 10 51 erg, M = IMq, and three progenitor 
. The + sign indicates the time at which 



radii, 10,3, 1 X 10 1 



Srn 



ph 



0.1. 



by C/O, eq. (TT31 is modified to 
T ph (i) = 1.5 eVf-° m7 2 



-0.2 p0.19 + - 
rt *,12 c 5 



0.35 



(27) 



In the absence of He, i.e. for Z = 1, T p h is simply given 
by eq. (p7]l. For a mixture of He-C/O, Z < 1, T ph may 
be obtained as follows. At high temperature, where He 
is still ionized, the He and C/O opacities are not very 
different and T p h obtained for a He envelope, eq. (|23l) . 
is similar to that obtained for a C/O envelope, eq. (|2T|) . 
At such temperatures, we may use eq. (|2"3")l for an enve- 
lope containing mostly He, and eq. ((27]) with Z = 1 for 
an envelope containing mostly C/O (a more accurate de- 
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Fig. 3. — Photospheric temperature (top panel) and the ratio 
of color to effective (=photospheric) temperatures (bottom panel) 
for explosions in a C/O envelopes (1:1 C: Q ra tio). The top panel 
compares the analytic approxi mat ion of cq. J27} (dashed lines) with 
the numerical solution of eq. (]_' I [ i using OP table opacities (solid 
lines). The color to effective temperature ratio was calculated as 
explained in § 13.21 Results are shown for E = 10 51 erg, M = 1 M®, 
and three progenitor radii, 10, 3, 1 X 10 11 cm. The + sign indicates 
the time at which 5m p ^ = 0.1. 




t, sec 



Fig. 4. — Photospheric temperature (top panel) and the ratio of 
color to effective (=photospheric) temperatures (bottom panel) for 
explosions in a He-C/O envelopes (1-Z=0.7 He mass fraction, 2:1 
C:0 ratio). The top pan el co mpa res the analytic approximations 
obtained usin g eq s. (1 23 P and 1271) with a transition temperature 
give n b y cq. 1281 1 (dashed lines) with the numerical solutions of 
cq. 1211 1 using OP table opacities (solid lines) . The col or to effective 
temperature ratio was calculated as explained in § 13.21 Results 
are shown for E = 10 51 erg, M = IMq, and three progenitor 
radii, 10, 3, 1 X 10 11 cm. The + sign indicates the time at which 
5m p h = 0.1. 



scription of the Z-dependence may be straightforwardly 
obtained by an interpolation between the two equations) . 
At lower temperature, the He recombines and the opac- 
ity is dominated by C/O. At these temperatures, T pn is 
given by eq. (l27l) with the appropriate value of Z. The 
transition temperature is given by 



r, 



He-C/O 



1Z 01 eV. 



(28) 



The photospheric radius, which is less sensitive to the 
opacity variations, is well approximated by eq. (I24[) . At 
the stage where the opacity is dominated by C/O, the 



bolometric luminosity is given by 



L = 4.7 x 10 



12 



p0.83 d0.8 
^51 n -*,l2 



/ p °- 14 Z°- 63 (M/M Q ) - 67 



£-0.07 ergs -l. 



( 29 ) 

The deviation of /„ from the assumed value of 2 for large 
6 m (discussed in § I2.2() does not affect significantly the 
results of eqs. (|2~7T) and (|28p. but may affect significantly 
the result given in eq. flM]). Eqs. (T7]l. (gUt and (27]) indi- 
cate that for an explosion with M = 10Af Q and E51 = 1, 
the luminosity will be reduced by a factor of ~ 2 (com- 
pared to the predictions of eq. (|29j) ) when T p h w 1 eV. 
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This further complicates the model for the emission on 
these time scales. See ^ 16. II for further discussion. 

In figures [3] and [4] we compare the analytic approxi- 
mation for T p h derived above to the results of numerical 
calculations using the OP opacity tables. For the C/O 
envelopes (figure [3]) , the approximation of eq. (|27|) holds 



to better than 6% down to T 



ph 



0.5 eV. For the Z=0.3 



mixed Hc-C/O envelopes (figure 2]), the approximations 
obtained by using eqs. (|23"|) and ([27]) with a transition 
temperature given by eq. ([25)) hold to better than 10% 
down to Tph — 0.8 eV. Using similar comparisons for dif- 
ferent compositions we find that similar accuracies are 
obtained over the range 0.7 > Z > 0.3, and for increas- 
ing or decreasing the C:0 ratio by an order of magnitude. 

3.2. Color vs. effective temperature 



We have shown in § 12.31 that photon diffusion is not 
expected to significantly affect the luminosity. Such dif- 
fusion may, however, modify the spectrum of the emitted 
radiation. We discuss below in some detail the expected 
modification of the spectrum. 

For the purpose of this discussion, it is useful to de- 
fine the "thermalization depth", rther> and the "diffusion 
depth", rdiff. r t hcr(i) < r ph(t) is defined as the radius 
at which photons that reach r p h(i) at t " thermalize" , 
i.e. the radius from which photons may reach the photo- 
sphere without being absorbed on the way. This radius 
may be estimated as the radius for which r sc tr a b s ~ 1 
(|Mihalas fc Mihalasl 11984) . where r sct and r a b s are the 
optical depths for scattering and absorption provided by 
plasma lying at r > rther(t)- ?"ther is thus approximately 
given by 

3(r t hcr - r ph ) 2 K sc t(rthor)Kab S (rthor)/3 2 (?"thcr) = 1, (30) 

where K sc t and K a b s are the scattering and absorption 
opacities respectively (typically, the opacity is dominated 
by electron scattering), rdiff is defined as the radius (be- 
low the photosphere) from which photons may escape 
(i.e. reach the photosphere) over a dynamical time (i.e. 
over t, the time scale for significant expansion). We ap- 
proximate rdiff by 



r-ph = r die 



+ Vci/ 3 '«sct(? , diff)p(? , diff), (31) 

where c is the speed of light. 

For rdiff < rthcr, photons of characteristic energy 
3T(r t honO > 3T p h will reach the photosphere, while for 
rthcr < ^diff photons of characteristic energy 3T(rdiff, i) > 
3T p h will reach the photosphere. Thus, the spectrum will 
be modified from a black body at T p h and its color tem- 
perature will be T co i > T p h- We approximate in what 
follows T co i = T(rther) for r d iff < rthcr and T co i = T(rdiff) 
for r d iff > rthcr- 

The lower panels of figures Q] [2] [3] and 0] present the 
ratio T co i/Tph for the various envelope compositions con- 
sidered. For this calculation, we have assumed that the 
scattering opacity is dominated by Thomson scattering 
of free electrons, used the electron density (as function 
of density and temperature) provided by the OP tables 
for determining K sc t, and estimated K a b s = k — n sct (re- 
call that k is the Rosseland mean of the opacity). It 
would have been more accurate to use an average of the 
absorptive opacities over the relevant wavebands, which 
are not provided by the OP table. However, since the 



dependence of the color temperature on the absorptive 
opacity is weak, T co \ cx , the corrections are not 

expected to be large. The figures imply that over the 
relevant time scale, t < 1 day, 



}t — T co i/Tph ~ 1.2 



(32) 



Using eq. (gSJ) with eqs. (T3J, ([23]) and ([27]) for the pho- 
tospheric (effective) temperature, the progenitor radius 
may be approximately inferred from the color tempera- 
ture by 



« 0.70 x 10 1 
for H envelopes, 



1 4 



Ccol 



(/t/1.2) eV 



.1.9 f 

z 5 Jp 



R* 



1.2 x 10 1 



T, 



col 



4.9 



1.9 rO.l 



j-1.9 J-U 

l 5 Jp 



(/t/1.2) eV 
for He envelopes with T > 1.07 eV, and 

5.3 



(33) 



(34) 



R* w 0.58 x 10 1 



T, 



col 



(/t/1.2) eV 



4- 9 / p U1 Z cm (35) 



for He-C/O envelopes when the C/O opacity dominates 
(the transition temperature is given in eq. (|28[l ). 

3.3. "Expansion opacity" 

We have neglected in our analysis the effective broad- 
ening of atomic lines due to the velocity gradients in the 
outflow. Line broadening may have a signi ficant effect 
on the opacity and on the observed emission flKarp et al.l 
119771: IWaeoner et "all 119911: lEastman fc Pintol I1993D . as 
wel l as on the dynamics (e .g. in the case of stellar winds, 
see lFriend fc Castorlll983T ). We give below a crude esti- 
mate of the line broadening effects for the problem of in- 
terest here. A detailed analysis of line broadening, which 
requir es detailed nu merical calculations (see, e.g., Sec. 
6.9 of lCastor|[200l . is beyond the scope of the present 
manuscript. 

The analysis of [Wagoner et all (|1991| ) shows that the 
effective line opacity introduced by the velocity gradi- 
ents may significantly affect the Rosseland mean opacity 
(at the relevant densities and velocity gradients, see their 
figs. 5, 6) at temperatures where recombination leads to a 
large reduction of the Thomson electron scattering opac- 
ity. The main contribution to this "expansion opacity" 
is from resonant line scattering of Fe group elements. As 
explained in § 13. 1[ the main effect that recombination 
introduces to our analysis is the penetration of the pho- 
tosphere to shells of high enough temperatures where 
significant ionization is maintained. Since the opacity 
enhancement due to velocity gradient effects does not 
prevent the strong reduction of the opacity due to re- 
combination, it will not prevent the penetration of the 
photosphere to a region of significant ionization. Nev- 
ertheless, the enhanced line opacity may introduce an 
order unity increase of the opacity at temperatures close 
to the recombination tempera ture at short, A < 0. 25 fi, 
wavelengths (see figs. 6, 7 of IWaeoner et al.l 119911 ). A 
detailed analysis of this effect is beyond the scope of the 
current manuscript. 

We should note, nevertheless, that the effective line 
opacity enhancement due to velocity gradients may be 
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smaller in our case compared to the estimates of ear- 
lier analyses. To show this, let us examine the following 
heuristic derivation of the effective broadening of atomic 
lines. Consider a photon that travels outward/inward 
in the region of relatively low optical depth. Due to the 
velocity gradient of the expanding ejecta, the photon fre- 
quency as measured in the plasma rest frame is shifted 
as it propagates by dv/v = —dv/c = —(dv/dr)dr/c = 
— (dv/dr)dt, where dv = (dv/dr)dr is the velocity dif- 
ference across dr and we neglect the plasma speed with 
respect to that of the photon in setting dr = cdt. Assum- 
ing that a photon is absorbed/scattered as its frequency 
is shifted across that of a line, then the probability for 
scattering/absorption is dP = dv(dN / dv) where dN/dv 
is the line "density" per unit frequency. The resulting 
photon mean free path is therefore l^ 1 ~ dP/cdt ~ 
c~ 1 \dv / dr\v{dN / dv). For v = r/t, which is valid at late 
time, we have 



i: 



\dv/dr 



■v(dN/dv) = {cty 1 v{dN/dv) (36) 



(compare, e.g., to eq. (3.10) of iWagoner et al.l [1991). 
The photon is absorbed/scattered provided the line op- 
tical depth is large enough. Neglecting the natural 
width of the lines, we may replace the line opacity with 
k v = kiVq5(v — vq). Denoting v' — (1 — dv/c)v we obtain 



Tl ~ / cdt'pK u > 



dv' 



v dv / dr 



dv /dr 



pui = ctpKi. 



(37) 



Thus, the line should be "counted" in determining dN/dv 
if ctpni ^> 1 (compar e, e.g., to eqs. (3 .10) and (2.7) of 
IWagoner et"atl (|1991h . and eq. (9) of iFriend fc Castorl 
(1983J) - 

The ana l yses of IWagoner et al.1 (|1991f l and of 
IKarp et al.l dl977l) are "local". That is, it is assumed 



there that as the photon's frequency is shifted by an 
amount comparable to the strong line separation, the 
parameters of t he plasma within which it propagate do 
not change (the lKarp et al.lfl977l assumptions are in fact 
more restrictive) . We have made the same assumption in 
deriving the final result of eq. (|37p . However, the valid- 
ity of this assumption is not obvious in our case. For the 
self-similar ejecta profiles described in § 12.11 p and T are 
steeply falling functions of Vf, p is roughly proportional 
and T is roughly proportional to vj 3 . Thus, as 



to v 



-in 



/ 



the photon moves outwards and its frequency is shifted 
(in the plasma frame) by ~ v/c (i.e. by dr/r ~ 1), p 
and T drop by factors of 10 3 and 10 1 respectively. This 
implies that r; may be significantly smaller than given 
by eq. §5$). 

4. REMOVING THE EFFECTS OF REDDENING 

We show in this section that the effects of reddening 
on the observed UV/O signal may be removed using the 
UV/O light curves. This is particularly important for in- 
ferring i?* , since i?* cx T" ol with 4 < a < 5 (see eas. [T3ll23l 
niul l27l: . 

The model specific intensity, fx, is given by 



T, 



aTp h g B b ( he/ AT col ) e " 



he 



(38) 



where 



gBB(x) = 



15 



1 



(39) 



D is the distance to the source, and t\ is the extinction 
optical depth at A. Let us define t\{t, A) by 

XT col [t = t x (t, A)] = A T co i(i), (40) 

for some chosen Ao- With this definition, the scaled light 
curves, 



/a [A, t\(t, A)] = 



D 



r ph(t\), 

x fx (A,i A ) 



T co i(t\) 



T ph (t x ) 







T co i(t\) 



(41) 



(where Tq is an arbitrary constant) are predicted to be 
the same for any A up to a factor e _TA , 

/a [A, t\(t, A)] = aT^g B B[hc/XoT col (t)} x e" T > . (42) 
he 

Let us consider now how the scalings defined above 
allow one to determine the relative extinction in cases 
where the model parameters {E, M, i?*} are unknown, 
and hence {T co i, T p h, r p h}(t), which define the scalings, 
are also unknown. For simplicity, let us first consider 
the case where the time dependence of the photospheric 
radius and temperature are well approximated by power- 
laws, 

rphoctf* T ph cx<- QT , (43) 

and the ratio Tcoi/Tph is independent of time. This is 
a good approximation for the time dependence of r p h in 
general, and for the time dependence of T co i and T p h for 
T ph > 1 eV (see eqs. [12] rjj [H M GEO and [32]). In this 
case eq. ([40]) gives 



tx{t,X)= 
and eq. (|41D may be written as 



1 / OiT 



t, 



(44) 



fx[\,t x (t, A)] = Const, x 
x/a A, 



A 
Ao 

l/a 

Ao 



(— 2a r +5aT) Jolt 



(45) 



The value of the constant that appears in eq. (|4"5)) . for 
which the normalization of fx is that given by eq. ([42]) . 
is not known, since it depends on the model parameters 
{E, M, i?*}. However, for any choice of the value of the 

constant, fx defined by eq. ([45]) is predicted by the model 
to be given by eq. (|42[) up to a wavelength independent 
multiplicative factor. Thus, the ratio of the scaled fluxes 
defined in eq. ([45]) determines the relative extinction, 



(46) 



/A[Al,tA(*,Ai)] 
A[A 2 ,tA(t,A2)] 



Let us consider next the case where the time depen- 
dence of T co \ and T p h is not a simple power-law. We 
have shown in § [5] and in § [3] that T co \ and T p h are de- 
termined by the composition and progenitor radius i?*, 
and are nearly independent of E and M. Adopting some 
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Fig. 5. — Ejecta velocity profiles at different times, pre- and 
po st- shock breakout, ob tained in the 14E1.3 model calc ulation 
of IBlin nikov ct al. (2000), compared with the [Matzncr & McKeel 
(1999) approximation (red line) for the post-breakout velocity pro- 
file used in our analytical model, given by eq. J3J for <5 m > 8 m go 
and assumed uniform at smaller <5 m (see eq. |6)l . 



value of eq. (l4*Tj)) may be solved for t\(t, X;R* 
eq. (I4ip may be written as 



and 



/a [A, t\(t, A; R*)] = Const, x r^T^tx^T^tx)- 1 

xA(A,t A ). (47) 

The model predicts therefore that scaling the observed 
flux densities using the correct value of i?*, the ob- 
served light curves at all wavelengths should be given 
by eq. ([4"2")) . up to a multiplicative wavelength indepen- 
dent constant. For this value of the ratio of the 
scaled fluxes at different wavelengths is independent of t 
and given by eq. (|46|). The value of i?* may be therefore 
determined by requiring the ratios of scaled fluxes to be 
time independent, and the relative extinction may then 
be inferred from eq. f4"o]l . We use this method in §[5] for 
determining i?* and the extinction curve for SN 2008D. 

5. COMPARISON TO OBSERVATIONS AND SIMULATIONS: 
RSG & BSC PROGENITORS 

5.1. SN1987A- A BSG progenitor 

Following the observations of SN 1987A, many numer- 
ical calculations modeling its light curve have been car- 
ried out (see e.g. lHauschildt fc Ensm an 1994). The lat- 
est and most comprehensive o f these calculations was 
carried out by IBlinnikov et al.l (|200Q[ ). and it provides 
UBV light curves from the time of bre akout to several 
month s following the explosion. The IBlinnikov et al.l 
(2000) radiation- hydrodynamics calculation, which in- 
cludes a detailed treatment of the opacities and a multi- 
group transport approximation for the propagation of 
radiation, should, to our understanding, capture all the 
relevant physics. 

In figures [5] and [Tj] we compare the results of our sim- 
ple m£dfiHo_those o f the detailed numerical calculations 
of IBlinnikov etaD (|2000h . We use the same progeni- 
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Fig. 6. — Comparison of early UBV measurements o f SN 1 987A 
with the 14E1 model calculations of Blinnikov ct al. (2000) and 
with our model calculations (red line) for similar model parameters 
(Eg_Y = 0.15 and a distance modulus of 18.5 assumed) . The 
upper most black curve in the figure (adapted from Blinnikov ct al. 
2000) is the V flux obtained in their 14E1A model calculation. 



tor parameters as those used in IBlinnikov et all ({2000) : 
A BSG (H envelope with n = 3, f p = 1) of radius 
R* = 3.37 x 10 12 cm, ejecta mass M = 14.67M Q , ex- 
plosion energy E = 1.03,1.34 x 10 51 erg for models 
14E1 and 14E1.3 respectively (and composition as in 
the outer part of th e progenitor given in figure 2 of 
IBlinnikov etHI [2000T> . Figure El compares the n umer - 
ical velocity profile with the iMatzner fe; McKeel ([1999D 
approximation we use in our model, given by eq. Q. 
The two agree to better than 10% over the relevant en- 
velope mass fraction. In figure [S] we compare the nu - 
merical early UBV light-curves of Bli nnikov et al.l (2000) 
with the ones calculat ed in our mo del, using eq. (|2H 
with the OP opacity (|Seatonl [2005). As can be seen 
in the figure, our analytic model gives fluxes which are 
larger by a factor of ~ 2 than those of the numerical cal- 
culation. This difference may be due t o differences be- 
tween the OP opacities and th ose used bvlBlinnikov et all 
(|2000f ). The opacity given in IBlinnikov et al.l (|1998D for 
p = 10- 13 gcnr 3 ,T = 15000°K and solar metalicity is 
larger than the OP opacity by roughly a factor of two. 
If a similar difference exists for the modified metalicity 
used in the SN1987A calculations, it would explain the 
luminosity discrepancy since the luminosity is roughly 
inversely proportional to k, see eq. (|15p . 



5.2. SNLS-04D2dc- A RSG progenitor 

The GALEX FUV and NUV measurements of the type 
Hp SNLS-04D2dc are shown in figure [7] Given the rela- 
tively low signal to noise ratio, we do not attempt here 
to constrain the progenitor parameters using a detailed 
analysis of the UV emission (as we do for SN2008D in 
§ [5]). Rather, we show that the observed UV flux is 
consistent with that expected from an expanding shock 
heated envelope of a RSG progenitor, and compare our 
simple model predictions to those obtained using detailed 
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Fig. 7.— Galex FUV (A cff = 1539A black) and NUV 
(^cff = 2316A magenta) observations of SNLS-04D2dc, not cor- 
rected for host and Galact ic extinc tion. The photometric anal- 
yses of IGezari et al.l J2008I) and of ISchawinski et al.l (120081 ) are 
shown in circles and squares respectively. The dashed (black, ma- 
genta) lines show the ( FUV, NUV) background levels inferred by 
Schawinsk i et al.l ( 12003 ). 

numerical calculations. For the latt er purpose, we use th e 
n umerical calculation s described in IGezari et al.1 ([2008). 

IGezari et al.l (|2008[ ) performed detailed numerical cal- 
culations, aimed at reproducing the early UV/O emis- 
sion from the SNLS-04D2dc, the progenitor of which is 
most likely a RSG. Their calculation was performed in 
two stages. At the first stage, a hydrodynamic calcu- 
lation of the explosion was performed using the one- 
temperatu re Lagrangian radia tion hydrodynamics code 
KEPLER (|Weaver et al.lll978l ). At the second stage, the 
emission of radiation at time t was calculated by solv- 
ing, usi ng the multi-group radia tion transport code CM- 
FGEN (jDessart fc Hillierl [2005h . the steady state radia- 
tion field for the hydrodynamic profiles obtained at time 
t, keeping the temperature profile fixed for r > 20 and 
allowing it to sclf-consistently change at smaller optical 
depths. 

In figure [8] we compare the IGezari et all 1)20081) cal- 
culations with the results obta ined by our m odel, using 
eq. dHJ) with the OP opacity (jSeatonl [20051 ) . The pro- 
genitor parameters we use are: A RSG (H envelope with 
n = 3/2 and f p = 25) of radius i?* = 865i? , explo- 
sion energy E$i = 1.44 and ejecta mass of M = 8.9M W . 
The st ellar radius is similar to that used by IGezari et al.l 
(2008), and f p = 25 was chosen (based on private com- 
munication) to provide an approximate description of the 
outer RSG envelope profile used in their calculation. The 
explosion energy and ejecta mass used in our model are 
20% larger and sm aller respectively than those used by 
IGezari et all (|2008f). i.e. the E/M ratio is 40% larger 
in our calculation. This value was chosen to reproduce 
the observed lumin osity. The light curves calculated by 
IGezari et al.l (pOOl) were shifted by ~ —1.5 magnitudes, 
i.e. the calculated luminosity was increased by a factor 
of ~ 4, to fit the observations. Since the luminosity is 
approximately proportional to E/M (see eq.H4|). this im- 
plies that for a given E/M ratio our calculation predicts a 




-2 2 4 6 8 10 

Fig. 8. — Comp arison of the numerical model calculations of 
IGezari et al. (2008) (thin lines) with our model calculations (over 
laid thick lines) for similar progenitor parameters and extinc- 
tion (see text for details). Due to the computationally demand- 
ing nature o f the n umerical calculation, the numerical model of 
IGezari et al. (2008) extends up to the time marked by the vertical 
dotted line (the thin curves at later times are a scaled SN II-P 
template). 

luminosity that is larger than that of lGezari et al.l (2008) 
by a factor o f ~ 3 (we have ve rified this by comparing 
the results of lGezari et a l. 2008, to our model results for 
the same E/M ratio). 

The following point should be noted here. Due to the 
large radius of the progenitor, the photosphere lies within 
the "breakout shell", i.e. 5 m , P h. < 5m,BO (see eq. H]), U P 
to t = tBo — 1.5 d (see eq. [H]). Our simple model is 
not valid at t < £bo- However, we expect it to provide a 
reasonable approximation for the photospheric tempera- 
ture and radius also at t < t^o, for the following reason. 
As long as the photosphere lies at 5 m , p h. < <5m,BOj the 
diffusion time at the photosphere is short compared to t. 
In this case we expect the spatial dependence of the radi- 
ation pressure to approximately follow par (see § 12 . 3[) , 
in which case the photospheric temperature is given by 
aTp h = 3p(r)/r. Since, as explained in § 12.31 the adia- 

batic pressure profile, p oc r , is similar to that obtained 
for short diffusion time, par, and since the adiabatic 
pressure profile is valid at all times for 5 m ^> 8 m ^o, 
we expect eq. (fl3|) to provide a good approximation for 
T p h also at t < ieo- The temporal dependence of r p h 
is expected to be somewhat stronger, at t < ieOj than 
the r p h oc t predicted by the simple adiabatic model, 
since the velocity profile at 5 m , p h. < ^ m ,BO is shallower 
than predicted by eq. . The largest deviation from the 
r p h oc < 8 behavior would be obtained assuming uniform 
velocity at & mtP h. < ^WbO; which would yield r p h oc t. 

The factor of ~ 2 discrepancy between the lumi- 
nosity predicted b y our model and that obtained by 
IGezari et al.1 (j2008) is due to the fact that our model 
predicts a somewhat, ~ 40%, larger velocity for the fast 
outer shells, and hence a larger photospheric radius. In 
a subsequent publication (|Rabinak et al.l l2010) we exam- 



12 





m 4 

I u 


CO 

1 


10 


E 




c_> 


-in 2 

10 


D3 


CO 


-in 1 
1 


7 


o 


10° 




> 

CD 


io- 








10" 2 




10 3 




Fig. 9. — A comparison of the temperature and density profiles 
(in blue and red respectively) of the ejecta given by the self-similar 
solution of §[2] (dashed lines) wit h those obtained in the "Harmonic- 
mean" model described in §[53] (solid lines), for an explosion of an 
n = 3 envelope with E51 = 6, M = 7Mq and = 10 11 cm, at 
t = 10 5 sec. 

ine the accuracy of the approximate ejecta density and 
velocity profiles described in § 12.11 for a wide range of 
progenitor models. 

6. THE EARLY UV/O EMISSION OF SN 2008D 

The analysis of the early UV/O emission of SN2008D 
is complicated by two major factors. Fir st, the extinction 
is larg e. It was loosely constrained by iSoderberg et al.l 
(2008) to 0.4 < E(B - V) < 0.8. This large extinction 
makes it difficult to extract the color temperature from 
the observations, and increases the uncertainty of the 
interpretation. The second complication arises from the 
fact that there is only one set of measurements in the UV 
at t < 2 days, and most of the relevant data points are at 
t > 1 d. This implies, as explained in § 13.11 that the sim- 
ple approximate solutions given by eqs. (|2"4"|) and (|2T|) for 
the photospheric radius and temperature are not accu- 
rate. For both He and mixed He-C/O compositions the 
reduction in opacity due to He recombination at > 1 d 
implies that the photosphere penetrates into the enve- 
lope beyond the range of validity of the approximation 
of eq. ([I]) for the initial density profile. For the analysis 
of SN 2008D, we will use therefore a more detailed de- 
scription of the density and pressure profiles of the ejecta. 
We first describe the modified ejecta profiles we use in 
§ 16. 1[ and then analyze the SN2008D observations using 
our model in § 16.21 As we show there, the deviation of 
the ejecta profiles from the S m —> self-similar profiles 
described in § 12. U lead to modifications of the model pre- 
dictions at t > 2 d. Comparison of our results to some 
earlier work appears in § 16.31 



6.1. Modified ejecta profiles 
The derivations of r p h and T p h in § 12.21 and in § [3] 



are based on the self-similar model of § 12.11 for the den- 
sity and pressure profiles of the ejecta. This model, in 
which p(vf) and p(vf) are both roughly proportional 



we adopt the " Harmon ic-Mean" model suggested by 
IMatzner fc McKeei (|1999fl . In this model, the density and 
pressure profiles are obtained by interpolating between 
the small 5 m (large Vf) self-similar power-law depen- 
dence of p(vf) and p(vf), and the power-law dependence 
p(vj) oc VJ, p(vf) oc Vf with a w 2 and f3 —1, ob- 
tained in the approximate analysis of Chevali er fe Sokerl 
(1989) for the lower velocity ejecta. This power-law de- 
pendence is obtained by assuming that the shock propa- 
gating within the ejecta may be approximately described, 
for large 5 m , by the Primakoff self-similar solution (a 
particular, analytic, case of the Sedov-von Neumann- 
Taylor solutions for shock propagation in to p oc r~" 
density profiles, obtaine d for uj = 17/7, e.g. !GafTetlll984t 
iBernstein fc Boo"k1 119801 ) . and by an approximate (self- 
similar) description of the post-breakout acceleration of 
the shocked plasma. The two power-law solutions de- 
scribing the larg e and small f f behay i or of t he density are 
matched in the IMatzner fe McKeei (|1999f ) "Harmonic- 
Mean" model at p = pbreak and Vf ~ u p brcak, and the 
pressure profiles at p = Phrsak and v f = ^nh rsak (see 
eqs. (46) and (47) of IMatzner fe McKeei [19991 ) . p hrc3 , k 
and fpbreak are determined by requiring the ejecta mass 
and (kinetic) energy to equal M a nd E respectively. 
jPhrsak and inbreak are determined in (Ma tzner fe McKeei 
1999) by examining numerical simulation results. We 
find that their parameter choice of pbrcak and inbreak 
leads an to over-estimate of the temperature, compared 
to that of the self-similar 8 m — > solution, by ~ 15% at 
S m = 10~ 3 . We therefore modify the value of pbrcak to 
obtain the correct self-similar behavior at 5 m — > 0. 

The " Harmonic- mean" density and temperature pro- 
files obtained as described in the preceding paragraph are 
compared in fig. [9] with those of the S m — > self-similar 
solution. Significant deviation from the self-similar pro- 
files is obtained for 5 m > 0.1. As discussed in cj 16.21 this 
deviation affects the model predictions for t > 2 d. The 
acc uracy of the " Harmonic- mean" model was examined 
in (jMatzner fc McKeei [T999) by comparing it to the re- 
sults of numerical calculations of the explosions of various 
RSG and BSG progenitors. Since it was found that this 
analytic model provides a good approximation for the 
envelope's profiles for different initial density structures 
of the progenitors, we expect the Harmonic-mean model 
to provide a good approximation also for the SN2008D 
envelope profiles. However, additional work, which is be- 
yond the scope of the current paper, is required in order 
to obtain a quantitative estimate of the accuracy of the 
approximation for S m > 0.1. 



6.2. Models 



observations 



to v 



-10 



is valid in the limit <L 



> (vf -> 00). 
In order to extend the model to larger values of 5 m , 



Spectroscopic observations have constrained the frac- 
tion of Hydrogen in the ejecta to < 5 x 10 -4 Mq 
(|Tanaka et al.l [20091) . We therefore consider below He 
and He-C/O envelopes. As explained in § [U the rel- 
ative extinction may be inferred from the light curves 
at different frequencies. However, for the clarity of the 
presentation we first analyze the data using two relative 
extinction curves that differ significantly in their short 
wavelength behavior, a Milky Way extinction curve with 
R v = 3.1 (here after MW) and a S mall Magellanic Clou d 
extinction curve (hereafter SMC) (|Cardelli et al.l[l~989l) . 
and only later show how the extinction curve may be 
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directly inferred from the data. 

Figure (|10[) presents a comparison of the color temper- 
ature T co \ and bolometric luminosity L inferred from the 
data with those obtained in our model for different pro- 
genitor and explosion parameters. We use the observa- 
tions of Swift/WTO (V, B, U, UVW1, UVM2, UVW2) 
Palom ar (g, r, i, z) (both t aken from iSoderberg et all 
12008ft . and FLWO (B, V, r, i) (IModiaz et al.ll2009ft . Since 
observations by different telescopes were carried out at 
different times, we interpolate the observations in differ- 
ent bands to times close to the observation times of the 
Swift /UVTO and the Palomar telescopes. In order to 
derive T co \(t) and L(t) from the observations, the mea- 
sured fluxes should be corrected for extinction. We carry 
out this correction by (i) assuming a specific extinction 
curve (MW or SMC) and (ii) determining the absolute 
value of the extinction by requiring T co \ inferred from 
the observations to agree with the model prediction at 
t ~ 2 d. This implies that for each set of model parame- 
ters (E, M, i?*, envelope composition), a different value 
of the absolute extinction, E(B — V) is chosen. 

For all the models shown in the figure we have used 
Er,i = 6 and M/M Q = 7, as suggested by Mazzali et al. 
(2008) from the spectral analysis of the observations at 
maximum light. As shown in § [2] and in § [31 E and 
M (mainly the ratio E/M) determine the normalization 
of the model luminosity, see eqs. (|25l) and ([29]) . but do 
not affect the time dependence of the luminosity, and 
T co \ is nearly independent of E and M, see eqs. (|23l) 
and (|27p . Since the conclusions drawn from the compar- 
isons in fig. (TTU)) are based on the time dependence of T co i 
and L, their validity is independent of the exact values 
chosen for E and M. 

Fig. (fTU]l shows T co i and L for models with three pro- 
genitor radii, i?* = 1,3,10 x 10 11 cm, and two com- 
positions, He dominated and mixed He-C/O composi- 
tion with a He mass fraction 1 — Z = 0.3. We show 
both the simple analytic approximations for L and T co \ 
given in § 13.11 (dashed curves) , which are based on the 
self-similar ejecta profiles of § 12.11 and L and T co \ ob- 
tained solving eqs. ([2"T]). (|50"|) and d3U) with the OP 
opac ity t ables for the modified ejecta profiles described 
in § 16.11 Since models with larger initial radii predict 
higher T co i, the absolute extinction inferred for models 
with larger radii is larger, E(B — V) = 0.625,0.7,0.8 
for i?* = 1, 3, 10 x 10 11 cm respectively. Once the abso- 
lute extinction is determined, from the comparison of ob- 
served and predicted T co i at t ~ 2d, T co \{t) and L{t) are 
inferred form the observations using the two relative ex- 
tinction curves (MW, SMC). In determining the observed 
T co \ and L, the photometric measurements are converted 
into monochromatic fluxes at the effective wavelength of 
the broad-band filters, and a BB temperature is deter- 
mined by a least-square fit to these fluxes. The resulting 
T co i and L are shown in fig. [TUJ The error bars represent 
the uncertainties obtained in the least-square fits. 

Examining fig. [10] we infer a small progenitor radius, 
i?* ~ 10 11 cm. Progenitors with larger radii require 
larger extinction to account for the observed flux dis- 
tribution at t = 2 day, which in turn imply that the ex- 
tinction corrected L(t) decreases with time at t < 2 d, in 
contrast with the roughly time independent L predicted 
by the models (see also eqs. ([25]) and ([29|)). and that the 



extinction corrected T co \ decreases faster than predicted 
by the models for t < 2 d. This is due to the fact that at 
earlier times the flux is dominated by shorter wavelength 
bands, for which the extinction correction is larger. Com- 
paring model predictions and observations at t > 2 d, we 
find that a mixed He-C/O composition is preferred over 
a He dominated one. The Z=0.7 model presented pro- 
vides a good fit to the data. We find that Z~ few 10's 
of percent is required to fit the observations. Note, how- 
ever, that the light curve at t > 2 d depends on the non 
self-similar part of the density and pressure profiles of 
the ejecta, for which we have used the " Harmonic-mean" 
approximation described in § 16.11 As mentioned in § 16.11 
additional work, which is beyond the scope of the current 
paper, is required in order to obtain a quantitative esti- 
mate of the accuracy of this approximation. We can not 
rule out, therefore, the possibility that the observations 
may be explained with a He dominated contribution and 
a density profile at large 5 m , that deviated from that 
given by the " Harmonic- mean" approximation. 

As explained in § [4] i?» and the relative extinction 
between different wavelengths may be determined from 
the O /UV light curves by requiring that the light curves 
observed at different wavelengths should all be given, af- 
ter scaling according to eqs. (fTO)) and (|3S"]l. by a sin- 
gle function, given by eq. (|42l) . In fig. [TT] we com- 
pare the measured specific intensities, corrected for ex- 
tinction and scaled according to eqs. (j40|) and (|45|) . 
with the model prediction, eq. (j42|) . For the scaling we 
have used {r p h(i), T co \(t), T p h(i)} obtained for the mixed 
He-C/O composition models presented in fig. [TU] with 
i?* = 10 11 cm and i?» = 3 x 10 11 cm. The extinction t\ 
was obtained by requiring the scaled intensity to best fit 
that predicted by eq. (|4"2")l (taking into account all data 
points at t < 4 d). The resulting extinction curves are 
shown in fig. (IT21) . As can been seen in the figure, the 
i?* = 10 11 cm model provides a much better description 
of the data than the R* = 3 x 10 11 cm model. The extinc- 
tion curve is more compatible with a MW extinction than 
with SMC extinction. It differs from the MW curve at 
short wavelengths, showing no prominent graphite bump. 
This, as well as the values obtained for Ay and E(B — V), 
Ay — 2.39 and E(B — V) = 0.63, are consist ent with the 
extinction inferred in ([Spdcrbcrg et al. 2008|, and private 
communication) . 

The following point should be explained here. We have 
obtained the absolute values of the extinction by adopt- 
ing some values for E and M. It is important to note, 
that, as explained in detail in §[H the relative extinction 
is independent of E and M, and may be inferred without 
making assumptions regarding their values. The model 
predicted temperature, T p h, and T co i/T p h are almost in- 
dependent of E and M, and depend only on i?* (and on 
the composition). E and M determine the normalization 
of r ph , r ph cx E 0A /M - 3 (see eq. ([24])). but do not affect 
its time dependence. Thus, modifying E and M changes 
the scaled fluxes of eq. |4T|) by some multiplicative factor, 
which is wavelength independent. Thus, the ratios of the 
scaled fluxes at different wavelengths are independent of 
E and M, and so are the inferred relative extinctions. 
Using the relative extinction inferred from the data, and 
adopting some relation between the relative and absolute 
extinctions, which determines the absolution extinction 



14 




> 

CD 



O 
CD 
CO 

CD 
i 

CD 



10 



10 



42 



10 



10 

t 



10 



I 






' — 3_ ' 







> 

CD 




10 
t 



10 



10 



10 

t 



10 



10 



44 



o 

CD 
C/5 

CO 
i_ 

CD 









- ^ 



> 

CD 



10" 



10 J 

t 



10° 




Fig. 10. — A comparison of the color temperature T co i and bolometric luminosity L inferred from the data with those obtained in our 
model for different progenitor and explosion parameters. Model results are shown for three progenitor radii, R t = 1, 3, 10 X 10 11 cm in the 
top, middle, and bottom panels respectively, and for two envelope compositions, He dominated (red) and mixed He-C/O composition with 
a He mass fraction 1 — Z = 0.3 (blue). E51 = 6 a nd M /Mq = 7 (and D = 27 Mpc) are assumed for all mode ls. S hown are both the simple 
analytic approximations fo r L a nd T co i g iven in ij l3.ll which are based on the self-similar ejecta profiles of § 12.11 (dashcdlines), and L and 
T co i obtained solving eqs. ( 1211 . J30D and H31I I with the OP opacity tables for the modified ejecta profiles described in § 16. II (full lines). The 
observed, extinction corrected L(t) and T co i(i) are inferred from the data by assuming a specific extinction curve (MW, green data points, 
or SMC, cyan data points), and determining the absolute value of the extinction by requiring T co i inferred from the observations to agree 
with model prediction at t ~ 2d. E(B — V) = 0.625,0.7,0.8 are inferred for ij* = 1,3,10 X 10 11 cm respectively. Blue data po ints are 
obtained for (best fit) extinction curves which are determined from the data itself, for the mixed He-C/O composition (see fig. 1121 1. 



(t\) we may therefore constrain the E 0A /M - 3 ratio by 
comparing the predicted and measured (absolute) flux 
at some frequency. Using the inferred E{B — V) = 0.63 
implies Ay = 2.39 for a MW extinction, for which we 
infer E 51 /(M/M Q ) = 0.75. 

6.3. Comparison with previous work 

Based on our analysis of the early UV/O emission 
of SN2008D we infer a small progenitor radius, i?» « 
10 11 cm, an E$i/{M /M©) « 0.8, a preference for a mixed 
He-C/O composition (with C/O mass fraction of 10's of 
percent), E(B — V) = 0.63 and an extinction curve given 
by fig- Q21 We compare below our conclusions to those 
of earlier analyses. 

Explosion models for SN 2008D were consid er by 
iMazzali et ail (f2008h an d b y iTanaka et all (f2009h . For 
E and M iMazzali et"all (|2008l) infer E ? i ~ 6 and 
M - 7M , i.e. E bl /{M/M Q ) ~ 0.85 while ITanaka et al.l 



(2009) infer E 51 = 6 ± 2.5 and M = 5.3 ± 1M Q , and 
E 51 /{M/M Q ) in the range 0.8 < E 51 /(M/M Q ) < 1.3. 
These values are consistent with our inferred value of 
E51 / (M/Mq) « 0.8. The progenitor radius inferred from 
the (jTanaka et al.l 120091) analysis is 0.9 < i?»/10 1:L cm < 
1.5, also consistent with our inferred i?* « 10 11 cm. 

In fig. [T3] we compare the photospheri c velocity ob- 
tained in our model with those obta ined bv lModjaz et al.1 
(2009) and by ITanaka et ahl pOOl analyzing SN2008D 
spectra. As can bee seen in the figure, our model predic- 
tions are in good agreement with the result s inferred from 
the sp ectral analyses. We also note that IMazzali et al.l 
(2008) infer, based on spectral analysis, that the mass of 
the ejecta shell moving at > 0.1c is ~ 0.03M©, consistent 
with our model prediction of ~ 0.02M© at this velocity. 

Let us discuss next our conclusion that the He en- 
velope contains a significant (10's of %) C/O fraction. 
This conclusion is consistent with the results of Maz- 
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Fig. 11. — The measured flux densities f\, corrected for extinction and scaled according to eq. 1141(1 (using D = 27 Mpc and To = 1 eV), 
compared with the model prediction of eq. (1421 (solid line) for the mixed He-C/O composition models presented in fig. llOl with R* = 10 11 cm 
and i?„ = 3 X 10 11 cm. Different symbols and colors show measurements at different bands: red circle— V, green x— B, blue left triangle— 
U, cyan right triangle- UVW1, magen ta square— UVM1, bl ack diamond- UVW2, magenta pentagon- r, red hexagon- g, black triangle- i, 
magenta right triangle- z, (taken from Sodcrbcrg ct al. 2008); black triangle- i', magenta diamond— r', red right triangle- B, green squarc- 
V (taken from lModiaz et aJjfeOOa) . 
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Fig. 12. — Comparis on o f the extinction curves inferred from the 
models described in fig llll (black and red points for R* = 10 11 cm 
and R* = 3 X 10 11 cm respectivel y) with Ay = 2.21 MW and 
A v = 2.16 SMC extinction curves (Cardelli et al.l[l989h . We find 



E(B -V) = 0.63,0.70 and A v 
Rt = 3 X 10 11 cm respectively. 



2.39,3.03 for R, = 10 11 cm and 



zali et al., who find a large fraction, ~ 10% of C at the 
fast v > 25000km/s He shells, rising to ~ 30% at the 
v ~ 20000km/s sh ells (P. Mazza l i, priv ate comm.). Both 
Mazzali et al. and lTanaka et all (|2009ft find a low, ~ 0.01 
mass fraction of O at the fast shells. 

Finally, we comment on the analysis of 
IChevalier fe Franssonl (|2008| ). who find a large pro- 
genitor radius, i?* ~ 10 12 cm based on both the 
X-ray and the UV/O emission. IChevalier fc Franssonl 
(2008) obtain T7* w 10 12 cm by interpreting the X-ray 
emission as thermal, Lx = 47r i? 2 erT^, and adopting 
Tx = 0.36 keV. Their motivation for a thermal interpre- 
tation of the X-ray emission, despite the fact that the 
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Fig. 13. — Comparison of the photosphcric velocity of the 
best fit model, in blue solid line, to results of other authors. In 
green squares are velocity interpreted from the He I A5876 line, 
and in green pentagon is the velocity from the analysis done via 
SYNO W to the "W" fea ture which later on vanishes both taken 
from[Modiaz ct al. (20(3). I n re d cir cles the photo s pheric velocity 
measured from spectral modeling by Tanaka ct al. (2009). 

spectrum is non thermal and extends beyond 10 keV, 
was that the non-thermal spectrum is difficult to explain 
theoretically. As explained in the introduction, the non- 
thermal spectrum is a natural conseq uence of the physic s 
of fast radiation mediated shocks (K atz et al.l [2010). 
Thus, there is no motivation, and it is inappropriate, to 
infer i?* assuming thermal X-ray emission. Moreover, it 
should be kept in mind that the X-ray breakout may take 
place within the wind surrounding the progenitor, at a 
radiu s significantly larger than (e.g. iWaxman et al.l 
2007). The discrepancy between their large radius and 
our small radius inferred from the UV/O data is due to 
the fact that we take into account the modification of 
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the opacity due to He recombination and the difference 
between color and effective temperatures. Since i?* oc n, 
neglecting the reduction of opacity due to HE recom- 
bination at T = 1 eV (t ~ 1 d) leads to a significant 
overestimate of the radius (compare eqs. ([551) and 
and neglecting the difference between color and effective 
temperatures leads to an additional over estimate of the 
radius, R* oc T 4 ol = (T col /T ph ) 4 T 4 h with (T col /T ph ) 4 » 2 
(see eq. (|34|l ). 

7. CONCLUSION 

We have presented a simple model for the early UV/O 
emission of core collapse supernovae. The photospheric 
radius, r p h(£), and (effective) temperature, T p h(i) are 
given for H envelopes by eqs. (fT2j) and ([T3|) in § and 
for He and mixed He-C/O envelopes (including pure C/O 
envelopes) by eqs. ([23]) -([28 ]) in § 13.11 T p h is determined 
by the composition and by the progenitor radius, i?*, and 
is nearly independent of the ejecta mass, M , and energy, 
E. M and E determine the normalization of r p h(i), but 
not its time dependence, r p h oc E 0A /M - 3 . The bolo- 
metric luminosity predicted by the model is nearly time 
independent (for T ph > 1 eV), see eqs. (14]), CLI]), ([25]) 
and (|29l) . Both r p h(i) and T p h(i) are only weakly de- 
pendent on n, the exponent determining the dependence 
of the progenitor's density on the distance from the edge 
of the star, see eq. ([T]). A discussion of the deviation 
of the spectrum from a black-body spectrum is given 
in § 13.21 where we find that the ratio of color to effec- 
tive temperature is approximately constant at early time, 
Tcoi /T ph « 1.2 (see figs. EE). 

For progenitor radii i?* < 10 12 cm, T p h approaches 
1 eV on day time scale. For He envelopes, significant 
recombination takes place at ~ 1 eV, leading to a sig- 
nificant reduction of the opacity, which, in turn, pre- 
vents Tph from dropping significantly below 1 eV, since 
the photosphere penetrates (deeper) into the envelope 
up to the point where the temperature is high enough 
to sustain significant ionization (see fig. [2]). Significant 
amounts of C/O in the envelope allow T p h to drop below 
~ 1 eV, since these atoms are partially ionized at lower 
temperatures as well (see fig. [J} . Model predictions de- 
pend only weakly on the C:0 ratio. 

Since T co \{t) is determined by the composition and by 
the progenitor radius may be inferred from the ob- 
served T co \. Eqs. ([33]) -([35" ]) give i?» as function of the 
observed T co \ for H, He, and He-C/O envelopes. A few 
comments should be made at this point. For a space and 
time independent opacity, i?„ oc kT 4 o1 (see eq. (fT^]) - ). In 
case k varies with temperature and density, the appro- 
priate value of k, i.e. its value at the photosphere at the 
time T C oi is measured, should be used. The fractional 
uncertainty in i?» is similar to the fractional uncertainty 
in k. The strong dependence of i?* on T co i implies that 
relatively small uncertainties in the determination of T co \ 
from the observations, or in its calculation in the model, 
lead to large uncertainties in the inferred -R*. Our ap- 
proximate treatment of the deviation from black-body 
spectrum due to photon diffusion implies T co \ is larger 
than T p h by w 20%. Estimating the uncertainty in the 
magnitude of this effect to be comparable to the magni- 
tude of the effect, implies a factor of w 2 uncertainty in 
the inferred i?» . A more accurate estimate would require 



a more detailed treatment of photon transport (including 
the effects of effective line opacity enhancement due to 
the large velocity gradients, see § !3.3p . 

Uncertainties in the observational determination of 
T co i are due to reddening. As explained in detail in § fU 
i?» and the relative extinction at different wavelengths 
may be inferred from the UV/O light curves. Scaling 
the observed light curves at different frequencies accord- 
ing to eq. (147)) , with T co i(t) and T p h(f) obtained in a 
model with the correct value of i?» , should bring all the 
light curves to coincide, up to a factor e~ TA where t\ is 
the extinction optical depth. The value of i?„ may there- 
fore be determined by requiring the ratios of scaled fluxes 
to be time independent, and the relative extinction be- 
tween two wavelengths may then be inferred from value 
of this ratio (see eq. (T4"6"10 . For the case where the time 
dependence of the photospheric radius and of the photo- 
spheric and color temperatures are well approximated by 
power-laws, which is a good approximation for the time 
dependence of r p h in general and for the time dependence 
of T co i and T p h for T p h > 1 eV, the relative extinction 
may be inferred independently of i?*, from the ratio of 
the fluxes scaled according to eq. (|4"5"]) . 

In § [5] we have compared our model predictions 
to observations of the early UV/O emission available 
for two SNe (SN1987A, SNLS-04D2dc), arising from 
RSG and BSC progenitors, and to detailed numer- 
ical (radiation transport-hydrodynamics) simulations, 
that were constructed to reproduce these observations 
([Blinnikov et all I2000t IGezari et all 12008ft . We have 
shown that our simple model may explain the observa- 
tions. We find, however, that our predicted luminosity 
is rs 2 times larger than obtained (for similar progeni- 
tor and explosion parameters) by the detailed numerical 
simulations. In the case of the BSG SN, this discrep- 
ancy is probably due to differences in the opacity tables 
we use and those used in the simulations (L oc 1/k, see 
§ 15.11) . In the case of the RSG SN, the discrepancy is due 
to the fact that our model predicts a somewhat (~ 40%) 
larger velocity for the fast outer shells, and hence a larger 
photospheric radius than that obtained in the numerical 
simulation. S ince t he details of the explosion model of 
IGezari et al.l (|2008[ ) are not given in their paper, it is 
difficult to determine the source of this discrep ancy. In 
a subsequent publication (jRabinak et al.l [2010ft . we will 
examine the accuracy of the approximate ejecta density 
and velocity profiles described in § 12.11 for a wide range 
of progenitor models. 

In § [5] we have used our model to analyze the early 
UV/O measurements of SN2008D. For this explosion 
we infer a small progenitor radius, i?» 10 11 cm, an 
E 51 /(M/M Q ) w 0.8, a preference for a mixed He-C/O 
composition (with C/O mass fraction of 10's of per- 
cent), E(B — V) = 0.63 and an extinction curve given 
b y fig . IT2l Our results are consistent (see discussion in 
§ I6.3[) with the E/M values inferred from modeling the 
light curve and spectra at maximum light ( Mazzal i et al.1 
2008J iTanaka et all 12009ft . with the i?* r ange obtained 
in the stellar evolution models described in lTanaka et al.l 
(2009), and with the extinction inferred from the analy- 
ses o f spectra at maximum light (e.g. ISoderberg et al.l 
2008). The photospheric velocity predicted by our 
model is consistent with the velocities inferred from spec- 
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tral analyses dMazzali et all 120081 iModiaz et all 120091: 
iTanaka et al.ll2009l ). see fig.fTl 

Our conclusion that the He envelope contains a sig- 
nificant (10's of %) C/O fraction is not as robust as 
the other conclusions, since it relies on observations at 
t > 2 d, where the emission is dominated by shells that 
were initially located at distances from the edge of the 
star for which the asymptotic (self-similar) description of 
the den sity, eq. (QJ , docs not hold (see detailed discussion 
in § 16. 2|) . In order to describe the light curve at t > 2 d we 
have used the approximation described in § 16.11 for the 
density and pressure profiles of the ejecta. Additional 
work, which is beyond the scope of the current paper, is 
required in order to obtain a quantitative estimate of the 
accuracy of this approximation. We can not rule out, 
therefore, the possibility that the observations may be 
explained with a He dominated composition and a den- 
sity profile that deviates from the approximation used. 
Nevertheless, our conclusion is consistent with the anal- 



ysis of Mazzali et al., who find a large fraction, ~ 10% of 
C at the fast v > 25000km/s He shells, rising to ~ 30% 
at the v ~ 20000km/s sh ells (P. Mazza l i, priv ate comm.). 
Both Mazzali et al. and ITanaka et al.l (|2009l ) find a low, 
~ 0.01, mass fraction of O at the fastest shells. 

The comparison of our analysis of the early emission of 
SN2008D with the analyses of the light curve and spec- 
tra at maximum light indicates that a combined model, 
describing both the early emission from the expanding 
and cooling envelope and the emission at maximum light, 
which is driven by radioactive decay, will provide much 
better constraints on the progenitor and explosion pa- 
rameters than those that may be obtained by analyzing 
either of the two separately. 



We thank P. Mazzali for useful discussions. This re- 
search was supported in part by ISF, AEC and Minerva 
grants. 
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